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Abstract 


A new method for design and generation of spiral bevel gears of uniform tooth depth 
with localized bearing contact and low level of transmission errors is considered. 

The main features of the proposed approach are as follows: 

(1) The localization of the bearing contact is achieved by the mismatch of the generating 
surfaces. The bearing contact may be provided in the longitudinal direction, or in the 
direction across the surface. 

(2) The low level of transmission errors is achieved due to application of nonlinear relations 
between the motions of the gear and the gear head-cutter. Such relations may be 
provided by application of a CNC machine. The generation of the pinion is based on 
application of linear relations between the motions of the tool and the pinion being 
generated. The relations described above permit a parabolic function of transmission 
errors to be obtained that is able to absorb almost linear functions caused by errors of 
gear alignment. 

A computer code has been written for the meshing and contact of the spiral bevel gears 
with the proposed geometry. The effect of misalignment on the proposed geometry has also 
been determined. Numerical examples for illustration of the proposed theory have been 
provided. 
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1 Introduction 


The research project is directed at the design and generation of face-milled spiral bevel gears 
with the following features: 

(1) The depth of the teeth is uniform which means that the height of the teeth is constant. 

(2) The gear tooth surfaces contact each other at every instant at a point. Thus, the 
bearing contact is localized and therefore the sensitivity of the gears to misalignment 
is reduced. 

(3) The surface contact under the load is spread over an elliptical area, whose center is the 
theoretical contact point. The set of instantaneous contact ellipses form the bearing 
contact. The developed methods of synthesis provide two types of bearing contact that 
may be directed: (i) in the longitudinal direction, or (ii) across the tooth surface. 

(4) Two generating surfaces, E f and £ c , are used for the generation of the pinion and 
the gear tooth surfaces and £ 2 , respectively. The dimensions of the instantaneous 
contact ellipse depend on the load applied to the gear drive, and on the chosen relation 
between the curvatures of the generating surfaces. 

(5) Gear misalignment may cause transmission errors of a high level and of such a shape 
that a high level of vibration will be resulted. The developed approach provides at each 
cycle of meshing a predesigned parabolic function that will absorb the transmission 
errors caused by misalignment. The cycle of meshing is determined as 


2tt 

91 = N~i 


(1) 
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where <j>\ and N\ are the angle of pinion rotation and the pinion tooth number. 


2 Method for Generation of Conjugated Pinion-Gear 
Tooth Surfaces 

Kinematic Relations 

We will consider initially an imaginary process for generation when the pinion and gear 
tooth surfaces will be generated simultaneously. Such an approach will permit one to obtain: 
(i) important kinematic relations to be executed on the CNC machine, and (ii) to visualize 

i 

the possibility to obtain two kinds of the localized bearing contact. In reality, the pinion 
and gear tooth surfaces are generated separately as it is discussed in the following sections. 

We start with the case when the axes of the pinion and the gear form an angle of 90° 
(fig. 1). However, the developed approach is applicable for gear drives whose shaft angle 
differs of 90°. The pinion and the gear during the process for generation perform rotation 
about axes X dl and respectively. Two rigidly connected generating surfaces, E, and E c , 
perform rotation about the Z m axis. Surfaces E f and E c generate the pinion and gear tooth 
surfaces £i and S 2 , respectively. The relation between the angular velocities w (1) and 
is represented by the equation 


ui 


u 


(‘) 

HI 


= sin 7i 


(2) 


where 71 is the pitch angle of the pinion. The transmission function of the gear drive 
must be obtained for each cycle of meshing as the sum of the theoretical linear function and 
the predesigned parabolic function as shown in fig. 2. We may represent the transmission 
function at the first cycle of meshing by the equation 
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( 3 ) 




Differentiating equation (3), we obtain that 


u/ 2 * = (— — 2a<pi)u>i 
■‘V 2 


Taking into account (see equation (2)) that 


(4) 


we obtain that 



sin7i 


(5) 


<t >2 = T J—r <i>c - a( / C ) 2 

iv2 sin 71 sin 71 


( 6 ) 


,(c) 

J*) = (ill - - 

7V2 sin 71 ' ' sin 7! 


, M 2a<j) c .. u 


(7) 


Considering 4> c and u c as the input parameters, we may obtain u>W, <f>\, <j> 2, and u/ 2) 
using equations (2), (5), (6), and (7), respectively. Since <f> 2 (<p c ) is a nonlinear function, 
the generation of the gear requires application of a CNC machine for the execution of the 


nonlinear function 4> 2 {4> c )- 

Note: The ratio of u; (1) and = u> (c) ) is constant during the process of pinion 

generation, the instantaneous axis of rotation is directed along the X m axis, and the pinion 
axode (in the process of meshing of the pinion tooth surface £1 with the generating surface 
£ ( ) is a circular cone with the pitch angle 


tan7i = 


Ni 

N 2 


( 8 ) 
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The axode of the generating surface is the plane Z m — 0 that is tangent to the pinion 
axode. 

( 2 ) ( 2 ) 

The ratio of angular velocities , and is not constant. Therefore, the axode of 
gear 2 being in meshing with the pinion and the generating surface £ c is not a circular cone. 
Generating Surfaces 

It was mentioned above that two rigidly connected generating surfaces £ t and £ c are 
applied for the generation of the pinion and gear, respectively. Two pairs of surfaces E t and 
E c are applied to provide two types of bearing contact for the generated pinion and gear 
tooth surfaces Ei and E 2 .. 

Case 1: 

The two generating surfaces are a cone £ t and a surface of revolution E c . The surfaces 
are in tangency along a circle and are rigidly connected each to other in the process for 
generation (fig. 3). Surfaces £ c and S 2 are in line tangency along lines L C2 while being in 
mesh in the process for generation. Similarly, £ t and £1 are in line contact along lines L t \ 
(not shown in fig. 3). However, instantaneous contact lines L C2 and L n do not coincide 
each with other but are in tangency at a point that belongs to the circle L ct . This means 
that surfaces £1 and £ 2 are in point contact at every instant that moves along the circular 
arc L ct in the process of meshing. The path of contact L ct (and the bearing contact) has a 
longitudinal direction. 

Case 2: 

The generating surfaces £ c and £ t are rigidly connected cones that are in tangency along 
their common generatrix L ct (fig. 4). Only contact lines L C2 between surfaces £ c and £ 2 are 
shown in the figure. Each contact line L C2 and the respective line L t \ are in tangency at the 
respective point of the generatrix L ct . This point is the current point of tangency of gear 
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tooth surfaces £ x and £ 2 . 

The generating surfaces E t and £ c are the surfaces of two head-cutters that are applied 
for the generation of the pinion and the gear, respectively. In reality, the head-cutters are 
provided with straight line blades or with circular arc edge blades but not with surfaces. Such 
blades are rotated about the head- cutter axes to form the generating surfaces. The angular 
velocity of rotation of blades must provide the required velocity of cutting or grinding, but 
is not related with the process for generation of the pinion or gear tooth surfaces. 

The head-cutter is installed on the cradle of the cutting machine and then performs with 
the cradle the rotation about the Z m axis with the angular velocity Details of 

the settings of the head-cutter on the cradle are discussed in sections 3 and 4. 


3 Generation of Gear Tooth Surfaces 

Applied Coordinate Systems and Machine-Tool Settings 

We have considered in Section 2 an imaginary process of generation of conjugate pinion- 
gear tooth surfaces based on the assumption that the surfaces will be generated simultane- 
ously. In reality, the pinion-gear tooth surfaces are generated separately, as it was mentioned 
above. 

This section deals with gear generation and the applied machine-tool settings. The gear 
tooth surfaces are generated by a head-cutter that is provided with straight-line blades (fig. 
5). Both sides of the gear tooth are generated simultaneously. The blades are rotated about 
the head-cutter axis. The edges of rotated blades form two circular cones as the generating 
surfaces used for the manufacturing of the gear. 

Henceforth, we will apply the movable coordinate systems S C2 and S 2 that are rigidly 



connected to the cradle and the gear, respectively (fig. 6). The fixed coordinate system S m 
is rigidly connected to the cutting machine. Coordinate system Sd 2 is an additional fixed one 
and is rigidly connected to the coordinate system S m . The orientation of Sd 2 with respect 
to S m is determined by angle 72 = 7 - 7i, where 71 is the pitch angle of the pinion pitch 
cone. During the process for generation, the cradle and the gear perform related rotations 
about the X m -axis and the A^-axis, respectively. The current angles of rotation of the 
cradle and the gear are designated by 4> C2 and <p 2 , that are related by equation (6) in which 
we change the designation of <f> c for <f> C2 . This relation provides the predesigned parabolic 
function designated for the absorption of transmission errors caused by misalignment. 

The installment of the head-cutter is shown in fig. 7(a). Figure 7(a) shows the position 
of the cradle when <p C2 =0, and the coordinate system coincides with S m . Coordinate 
system S n is an additional coordinate system that is rigidly connected to S CJ , and performs 
rotation with S C2 as shown in fig. 7(b). Coordinate axes of 5^ and S C2 have the same 
orientation. Axis Z V2 is the axis of the head-cutter (see below). The installment of the 
head-cutter in S C2 may be determined by the parameters of the triangle with the apexes 
0 C2 , 0 P2 , and M, where M is the mean point of the cone distance of gear 2 (fig. 8). The 
installment of the head- cutter in coordinate system S C2 may be determined as well by Hg 
and V G (fig- 8), where 

H a = A m — FL 2 sin 0 2 (9) 

V G = R u 2COS& (10) 

Here, R u2 = \0 P2 M\ is the nominal radius of the head-cutter (fig. 8), A m = \0 C2 M\, 
and /?2 is the spiral angle of the gear. Alternative coupled parameters of installment, S r2 = 
| Oc 0^ | and q 2 , Hg and Vq, are related by the equations 
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( 11 ) 


St 2 = {Hi + V^) 0 - 5 

92 = tan- 1 (-^) (12) 

Equations of Gear Generating Surfaces 

We have mentioned before that two cones as generating surfaces are used for cutting of 
the space of the gear (fig. 5). Fig. 9(b) shows the cutting blade, and fig. 9(a) shows one 
of the generating cones. The generating cone is formed by rotation of the blade edge about 
the Z V2 axis. The equations of the cone are represented in S V2 as 


T n {S a ,0 G ) 


(r c — Sg sin og) cos Qq 
(r c — Sg sin arc) sin 6g 
—Sg cos c*g 


( 13 ) 


where 0 G and S G (fig- 9(a)) are the surface coordinates (the Gausian coordinates); a G is 
the blade angle; r c is the radius of the head-cutter that is measured at the tip of the blade. 


Parameters r c and /? u2 are related as follows 


r c = Ru 2 ± 


PW 

~2~ 


(14) 


The upper and lower signs in (14) correspond to the gear concave and convex sides. 
Equations (13) may represent both generating cones, if we will use the rule of signs in 
equation (14), and consider that a G > 0 and a G < 0 correspond to the gear concave and 
convex sides, respectively. 

The unit normal to the gear generating surface is represented by the equations 


„ =T£3_ N = 9rn :: dTf * 

” |N P2 | ’ " Me dSc 

Equations (14) and (15) yield 


( 15 ) 
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(16) 


n n (^) 


— cos cxq cos 9q 

— cos og sin Oq 

sin etc 


Algorithm for Derivation of Gear Tooth Surfaces 

The derivation of gear tooth surface S 2 is based on the following procedure: 

Step 1: 

Initially we derive the family of generating surfaces in S 2 using the matrix equation 


r 2 {S g ,0g,4>c7) = M 2 p 2 rp 2 (5 G , 0 G ) (i 7 \ 

= M 2 M 2 )M dim M mC2 (<j> CJ )M cjn r P2 (S G ,e G ) K ’ 

where <j > 2 and 4> C2 are related by equation (6). 

Step 2: 

Then, we derive the equation of meshing between the gear and the generating surface 
that we represent in the form 


nl? ) -v i * 2) = f(S G ,0G,<i>ci) = 0 (18) 

where n^ 2 * is the unit normal to the generating surface S P2 , and v (pi2) is the relative (sliding) 
velocity. The scaler product (n^ P2 ^ • v^ P22 ^) may be represented in any coordinate system. In 
our derivations, the scaler product is represented in S m . 

Note: we may represent as well the equation of meshing in the form 

N£ 2) • v( p22 > = 0 (19) 

where is the normal but not the unit normal to the generating surface. 

Equations (17) and (18), considered simultaneously, represent the gear tooth surface as 
the envelope to the family (17) of generating surfaces. 
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Derivation of Family of Surfaces (17) 


Vector function t P2 (Sg,8g) has been already represented by (13). Matrices in equation 
(17) are represented as follows (see figs. 6,7, and 8) 


M 2 d, = 




1 

0 


0 


0 0 


MtoC2 — 


0 

— sin ' 

0 

COS 4> C2 

sin <t> C2 
0 
0 


1 

0 

0 


sin <p 2 

0 

cos <f> 2 

0 

0 

1 . 

sin 72 

0 ' 

0 

0 

COS 72 

0 

0 

1 . 

1<t>C2 0 

0 ■ 

<Pc 2 0 

0 

) 1 

0 

) 0 

1 


( 20 ) 


( 21 ) 


(22) 


Equations (17) and (20)-(22) represent the family of generating surfaces in form of three 
parameters. The relation between these parameters is represented by the equation of meshing 
(18). 

Derivation of Equation of Meshing 

We represent in coordinate system S m the unit normal nj£ 2 * to the generating surface 
using the following matrix equation 


n g>>(»e,*) = I'mC2 ^* C 2P2 ^P2 (»c) 


(23) 


Here: 

Vector function n P2 (0c) is represented by equation (16). Matrix L is the 3x3 submatrix 
of the 4x4 matrix M. After derivations, we obtain 


rim(^G? 4*02 ) 


— COS Qg COs(0G + <t>c 2 ) 

- cos Qg sin(#G + <j>c 2 ) 

sinQc 


(24) 


11 



As a reminder, the variables fa and ^ are related by equation (6). 

The derivation of the relative velocity is based on the following considerations 


Step 1: 

We represent vector v^ 22 ^ as follows 

vS? 2) = (w<P» - x r<?> (25) 

where 

u>£ 2) = [0 0 cj (ci) ] t ( 26) 

= L m2 u> (2) = L md2 L <i2 2W (2) = [u» (2) cos 72 0 u> (2) sin7 2 ] T (27) 

The relation between w (p2) and u; (2) may be obtained by differentiating equation (6) that 
yields 


JV _ w lc 2 ) _ 2 a { 4 ^-)^- (28) 

N 2 sm 71 sin 71 sin 71 

4 Machine-Tool Settings for Pinion Generation 

Introduction 

Henceforth, we will consider methods for the pinion generation that provide two types of 
the path of contact : across the surface and along the surface. The developed machine-tool 
settings will be determined: (i) for the imaginary process for generation based on application 
of two rigidly connected generating surfaces (see figs. 3 and 4), and (ii) machine-tool settings 
represented in terms of the Gleason Works terminology. 
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Coordinate Systems Applied for Pinion Generation 

Movable coordinate systems S C1 and S\ are rigidly connected to the cradle of the cutting 
machine and the pinion, respectively (fig. 10). An auxiliary coordinate system S Pl is rigidly 
connected to the head-cutter and the cradle. Axes of coordinate systems S Cl and S Pl have 
the same orientation. Fixed coordinate systems S m 

and Sdi are rigidly connected to the cutting machine. The pinion and the cradle perform 
rotation about axes x dl and z m , respectively. During the generation, the cradle and the 
pinion are rotated uniformly, and the angles of rotation are related as follows : 

<t> ci = <£1 sin 7! (29) 


where 71 is the pinion pitch angle. 

Pinion Machine-Tool Settings 

The pinion machine-tool settings for the imaginary process of generation are determined 
by the following set of parameters (fig. 10): 71, S r \ and qi- The machine-tool settings for the 
real process of generation may be determined by turning of systems 5 1 , S dl , S Pl and S Cl at 
180° about the x m -axis. Fig. 11 shows the installment of the pinion on the cutting machine 
and its rotation during the process for generation. Fig. 12(a) shows the initial installment of 
the head-cutter and cradle and fig. 12(b) shows the rotation of the cradle during the process 
for generation (fig. 12(b)). 


5 Equations of Pinion Tooth Surfaces 


Orientation of Bearing Contact Across Surface 


13 



The mentioned type of bearing contact is provided by application of two rigidly connected 
circular cones that generate the gear and pinion tooth surface, respectively (fig. 4). Fig. 13 
shows the cross sections of the surfaces of the pinion and gear head-cutters in plane z m — 0 
when the cradle is at the initial position that is determined with <f> Cl = 0, and the coordinate 
systems S C1 and S m coincide each with other (fig. 10(b)). 

Fig. 14 shows the profiles of the cutter blades used for the generation of the pinion 
concave side and convex side surfaces, respectively. Fig. 15 shows the generating surface 
(cone) in system S P1 . The derivation of the generated pinion tooth surface is based on the 
following procedure. 

Step 1: 

The head-cutter generating surface is represented in S P1 as 

(jcp — S p sin Op) cos 6 P 
(rep - S p sin q p ) sin 0 V 
—S v cosa p 
1 

Here, is the point radius (fig. 14), a p is the blade angle. The sign of a p should be 
considered as positive and negative when the pinion convex and concave tooth surfaces are 
generated, respectively. Parameters 0 V and S p are the head-cutter surface parameters. 

The unit normal of the head-cutter generating surface is represented in S Pl as 





(31) 


( 32 ) 


14 



The family of head-cutter generating surfaces is represented in S\ by the matrix equation 


r l(*$p> Qpifici ) — MidjMrfjmMrncjMcjpj r pi (i9p, 0p) 


(33) 


Here, 


Mh, = 


10 0 0 
0 cos <f>i sin 0 
0 — sin <f> i cos <f>\ 0 


^1 d ] m — 


Mmci 


o 

o 
1 

0 

1 J 

cos 71 

0 sin 7 a 

0 1 

0 

1 0 

0 

— sin 71 

0 cos 7 i 

0 

0 

0 0 

1 . 

’ cos<j> Cl 

— sin <t> Cl 0 

0 

sin <p Cl 

cos <f> Cl 0 

0 

0 

0 1 

0 

0 

0 0 

1 

ri 0 

0 Sri COS qi 

1 


M = 

iViciPl 


0 10 5 r isin9i 


0 0 1 
0 0 0 


0 

1 


Step 3: 

The equation of meshing is represented in system S m as 


(34) 


(35) 


(36) 


(37) 


n • v (cil) = 0 

ll m v m 


(38) 


Here, 


n m {d p ,<? c: ) = 


rm 


,l ym 

^:m 


— Lmcx ^cipi ^ci (^p ) 


V^ l}) = x r„ 


(39) 

(40) 


where 
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( 41 ) 


— Mmci ^^cipi ^pi 


5 P sin a p cos (6 P + <t> Cl ) + r cp cos(0 p + <j> Cl ) + S T i cos(<£ Cl + 91 ) 
— S p sin a p sin(0 p + <£ Cl ) + r cp sin(0 p + <j > Cl ) + S r i sin(<^ Cl 4- 91) 


—S p cos a p 


Using the designations 


Bi = r cp cos(0 p 4- <f > Cl ) + ^ri cos(<£ Cl + 91) 1 

f? 2 = rcp sin(0 p + <i > ci ) + S r i sin(<£ c , + qi ) J 


(42) 


we obtain 


fm — 


—S p sin a p cos (0 P + <t>a ) + B\ 
—S p sin a p sin(0 p + <£ C i ) + B 2 


— S p cos a p 


(43) 


w; 


(cil) _ w («i) _ ^,(1) = 


' 0 


’ U>1 COS 7 1 


UJi cos 7 j 

0 

— 

0 

= 

0 

- . 


u?i sin 71 


0 _ 


Then, we obtain that (since u ci = u>\ sin 71) 


(44) 


0 

v^ 11 ' = u;^ 1 * x r m = S p cos a p cos 71 

_ -S p sin a p sin(0 p + <f> Cl ) cos 7! + B 2 cos 71 J 

We have assumed not loosing the generality of the approach that |u>i| = 1. 
The equation of meshing yields 


(45) 


n ym S p cos q p cos 7i — n xm S p sin a v sin(0 p + <f > Cl ) cos 71 + n zm B 2 cos 71 — 0 (46) 


or 


Sp{0 p , ) 


-n zm B 2 

n ym cos Op - n zm sin a p sin (6 p + o cl ) 


(47) 


Step 4: 
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The envelope to the family of generating surfaces may be represented in two parameter 
form if we eliminate in equation (33) parameter S p using equation (47). 

Orientation of Bearing Contact in Longitudinal Direction 

The mentioned orientation of the bearing contact can be achieved by application of 
generating surfaces shown in fig. 3. The profiles of the blades of the head-cutter used for 
the generation of the pinion are shown in fig. 16. 

The derivation of the equations of the pinion tooth surface generated by a head-cutter 
with circular arc blades is based on the following procedure: 

Step 1: 

The coordinates of the center of the blade circular arc for the concave side are represented 
by the following equations (fig. 16) : 


0 0 C X = OM x + M X C X 


(48) 


where 0M X = r Cl , M X C X — R\ 

Equation (48) yields the following equations for the coordinates of center C x 


x ( a c,) = 0 0 C X • i 0 = r ci - R x cosQp 

y£ Cl) = OpCl Jo = 0 

z (Ci) - Q o C 1 . lc 0 = -R x sinap 

Step 2: 


(49) 


The position vector of a current point of the circular arc is represented in S 0 by the 


equation (fig.l 6) 


0 o A = 0 0 C X + C X A = 0 0 C X + R x n 0 


(50) 
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where 


n 0 = [cosAi 0 — sinAi] T (51) 

is the unit normal to the circular arc that is represented in S 0 . 

Step 3: 

The head-cutter surface is a surface of revolution that is generated while the circular 
arc is rotated about the axis of the head cutter. The head-cutter surface is represented in 
coordinate system S V1 as follows 


r p , = r'f‘1 + R t n n 


(52) 


Here 


_(Ci) = L r (Ci) 
F P1 lj PlO l o 


is the position vector of point C\ that is represented in S\ 


npj — Lpjoiio 


is the unit normal to the generating surface that is represented in S VI 


L P io 


cos d p — sin dp 0 
sin dp cos dp 0 

0 0 1 


(53) 


(54) 


(55) 


Step 4: 

For the following derivations, we represent the surface of the head-cutter in coordinate 
system S m (fig. 10(b)). We may use for this purpose the following vector equation 


r m = 0 C1 0 P1 + + #in r 


(56) 
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Here (fig. 10(b)) 


0 Cl Op, = S T 1 [cos qi sin qr 

0] T 

(57) 

r m l} = LmC^p, 


(58) 

n m — h m c, rip. 


(59) 

cos 4> c — sin <f> c 

0 ' 


LmC, = sin <t>c cos <j> c 

0 

(60) 

1 

O 

o 

1 



Step 5: 

The pinion tooth surface is represented in as the envelope to the family of the head- 
cutter surfaces that is generated in coordinate system S\. The equation of the family of 
surfaces is as follows : 

ri(Ai,0 p , 4 > Cl ) = M 1 i(^i)M (fim r m (Ai,^p,^ Cl ) (61) 

As a reminder, angles <j> Cl and <t> i , angles of rotation of the cradle and the pinion, are 
related by the equation 

<f>a = <t> 1 sin 7i ( 62 ) 

Step 6: 

The envelope to the family of surfaces is determined by equation (61) and the equation 
of meshing that we represent as follows 

■ n m = 0 (63) 

The final expression of the equation of meshing is based on the following derivations. 
The sliding velocity is determined as 
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v£ pi) = w£ pi) X r m = (w£> - «<»>) X M mpi r Pl = u£ Pl) x (r£ lJ 4- fl,n m ) (64) 
Equations (63) and (64) yield 

(w£ w) x r£ l} ) ■ n m = v£»^> (65) 

where 

tf' 1 = L mpi r pi (66) 

and vJ* pl ’ Cl ) designates the relative velocity for point C\. 

Step 7: 

The surface unit normal n m is a vector function of three variables (Ai,0,<£ c ). We may 
simplify the vector function n m (Ai,#, <j> c ) using the following considerations : 

(1) Equations 


v^>n m =0 (67) 

and 

vS " c '' -».=0 ( 68 ) 

yield that the relative velocities determined at the point of tangency of the head-cutter and 
the pinion, and at the center of the circular arc are collinear. Taking this into account, we 
may represent the unit normal n m to the surface of the head-cutter by the following equation 


X y(lPl.Cl) 


\Tm X V m 


(lpi £\) | 


(69) 
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where T m is a unit vector of the tangent to the 0 p coordinate line on the head- cutter 
surface, that is represented in S m system. The unit vector T m is represented by the following 
equation 

= L mpi L pi0 j 0 Lmpi^pi 

Vector j 0 is perpendicular to the plane ( x 0 ,z o ) of the circular arc (fig. 16). A point of 

the circular arc traces out in S P1 a circle, and the unit tangent t pi is represented as 

— sin 0 P 

T Pi = L Pl0 j 0 = cosOp (71) 

0 

The advantage of application of equation (69) is that we may represent the surface unit 
normal by the vector function n m (0 p , <f> Cl ). Then, the equation of meshing will yield the 
relation 

v£ p, ’ Cl) ■ n m = f(0 P ,4> Cl ) = 0 (72) 

that is free of parameter Aj. 

Angle Aj can be obtained from the equation 

cosAi = r m • n m (73) 

Note : Similarly, we may derive the equations of the convex side of the pinion tooth 
surface. Center C 2 of the circular arc (fig. 16) is represented in S 0 by the equations 

x i C2) = r C2 + R-2 COS Ctp 

=0 [ (74) 

z (Ci) = -R 2 sina p ; 
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6 Computerized Simulation of Meshing and Contact 

The goals of this investigation are : (i) the determination of paths of contact for aligned 
and misaligned gear drives, (ii) and the determination of influence of misalignment on the 
transmission errors and the shift of the bearing contact. 

Conditions of Continuous Tangency 

We set up three coordinate systems S fc , S x , S 2 (fig. 17) that are rigidly connected to 
the frame, the pinion and the gear; <j>i' and <j> 2 are the angles of rotation of the pinion 
and the gear when they are in mesh; H, V, Q 8' are the parameters used to simulate the 
misalignments that represent. 

The contact of the tooth surfaces is localized and they are in tangency at every instant 
at a point. The simulation of meshing is based on the condition of continuous tangency of 
pinion-gear tooth surfaces £ x and E 2 , that are represented in coordinate system Sk as follows 

= ( 75 > 

= ( 76 ) 

where |n^| = |n^ | = 1; (9 P , <f> ei ),($c,<j> C2 ) are the surface parameters of the pinion and 
the gear, respectively; <j>\ <f> 2 are the angles of rotation of the pinion and the gear being in 
mesh. 

Equations (75) and (76) represent a system of five nonlinear equations in six unknowns 
represented as 


fi{0p,<f> cl ,<p\,0G,<f>c2^2) ~ 0 (* = 1 > 5 ) 


(77) 


where /, € C' 
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One of the unknowns in equation system (77), say <f>\, is chosen as the input one. The 
continuous solution of these equations is an iterative process that is based on the following 
procedure. 

Using the first guess, we consider that a set of parameters designated as 

p(°> = («w, (£)«■>, (£)'•>) (78) 

satisfies equation system (77). We assume as well, that we have 

A s = + o (79) 

Then, in accordance to the Theorem of Implicit Function System Existence [15], equation 
system (77) can be solved in the neighborhood of P ^ by functions 

&M), K Wi ). Mrfi). K (#), 4>'M\) (8°) 

The solution of these nonlinear simultaneous equations is found by numerical methods. 

The path of contact on the pinion tooth surface may be represented by the following 
expressions 


ri(0p,^c), 9p(4> i). <i>cM i) 

Similarly, the path of contact on surface S 2 may be represented by 


(81) 


r 2 (Qg,4>c 3 ), e G m, i) ( 82 ) 

The transmission errors caused by misalignment are determined by the equation 


A </>2 — d 2 N 2 < ^ 1 


(83) 
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Derivation of Equation System (80) 

We use for derivation matrix equations (84)— (87), equations of surfaces E v and Eg* and 
the surface unit normals. 


4 l) (»* 


= M«W)r 

1 0 0 
0 cos <t>\ sin <f>\ 
0 — sin <f>[ cos <p\ 
0 0 0 


H 

0 

0 

1 


ri W P ,K) 


(84) 


d 2) ( 0 c, V*, <? 2 ) = M w (^)r 2 {0c, 4 >c 2 ) 


cos (7 + 5') 0 - 

- sin (7 -1- S') 

Q 1 


1 0 

0 

0 ‘ 

0 1 


0 

-V 


0 cos <j> 2 

— sin <p 2 

0 

sin (7 4- o') 0 

cos 

(7 + V) 

0 


0 sin <f >2 

cos <j> 2 

0 

0 0 


0 

1 


0 0 

0 

1 . 


' 1 

0 

0 






0 

cos 4>\ 

sin o' 

t 

1 

ni (0 P ,<t> cl ) 




0 

— sin <f>[ 

cos 6 

/ 

1 . 





r 2 {0a, 4>a) 


(85) 


( 86 ) 


n 1 2 ) ( 0 g>c 2 , o' 2 ) 



cos (7 + S') 0 — sin (7 + 6 ') 


'10 0 

— 

0 1 0 


0 cos <£2 — sin q >2 


sin (7 + S') 0 cos (7 + S') 


_ 0 sin (f>' 2 cos <^>2 


n 2 (0g,4*) 
(87) 


Here; 71 is the pitch angle of the pinion, 72 = 7 — 7i> 7 is the angle formed by the 
pinion-gear axes of rotation. Usually, 7 = 90°. 

Be advised that in the case when the generating surface is a cone, the unit normal is 
represented by a vector function of one variable (see Section 3) 


Bearing Contact 

Theoretically, the pinion and gear tooth surfaces are in point contact. Under the load, 
the contact is spread over an elliptical area. The determination of dimensions and orientation 
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of the instantaneous contact requires the knowledge of principal curvatures and directions 
of the contacting surfaces. The solution to this problem is substantially simplified due to 
representation of the curvatures of the generated surface by the curvatures of the generating 
surface and the parameters of motion [15]. 

For the case when the pinion generating surface is a cone, the principal curvatures of the 
pinion generating surface are represented as 

= 0 , ( 88 ) 

= cos Op/ (r cp - S p sin a p ) 

The principal directions on the pinion generating cone are 


ij 1 ) = [— sin 6 P cos 6 P Op 


.(») - 


= [sin q p cos 9 p sin a p sin 0 p — cos a p ] 


(89) 


For the case when the pinion generating surface is a surface of revolution, the principal 
curvatures of the pinion generating surface are represented as 


k ( ' ] = 1 /Rx 

k { 2 1] = cos + Ri ■ cosAj ) 

The principal directions on the pinion generating cone are 


(90) 


= [— sin B p cos 0 P 0] T (g-^ 

= [sin Ax cos # p sinAiSin^p — cosAi] r 

Similarly, when the generating surface of the gear is also a cone, we have that the principal 
curvatures of the gear generating surface are 


Jtj 2) = 0 

= cos otQj(r c — Sg sin ag) 
The principal directions on the gear generating surface are 


(92) 
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( 93 ) 


ef* = [— sin Og cos Og 0] t 

= [sin qg cos Og sin q<j sin Og — cos qtg] T 

The principal directions are represented in systems S Pl for the pinion case, and S n for 
the gear case, respectively. 

Because the generating and generated surfaces are in line contact, their principal curva- 
tures and directions are related by the following three equations [15] 


tan2cr 


2£ 13^23 

^23 ^13 (^/ kh)t 33 


kn — k, = 


- 2 * l&Z 
f 33 sin 2 a 


(94) 


k q + k, = kf + h+ ^t ^ 

l 33 

Here, kj,k^ are the principal curvatures of the generating surface; k s , k q are the principal 
curvatures of the generated surface; cr is the angle between the principal directions of the 
generating and generated surfaces (fig. 18(a)). The expression for i i 3 , ^ 23 i and t 33 are 
represented in [15]. 

The determination of the dimensions of the contact ellipse and its orientation is based 
on application of following equations 


where 


tan2r = 


g 2 sin a 


( 12 ) 


9 \ -92 cos a 


( 12 ) 


(95) 


91 = k[ l) - k[ 2) 

92 = k^ ] - k\ 2} 
a (!2) _ cos - 1 (e^ • ef^) 
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j fcj 1 * and are the principal curvatures of the pinion surface. k\ 2 ’ and k^ are the 

principal curvatures of the gear surface. 

The major axis and minor axis of the contact ellipse may be determined as 

2«=2 26=2 ( 96 ) 

where S is the elastic approach obtained from experimental data; A and B are determined 
by 


A = “ k x * “ (Si “ 2 0 i 0 2 cos 2a(12) + 02^1 ('97') 

B = - k§ ] + (5? - 2^2 cos 2 a< 12) + <?|)2] 

and 

di) _ d i) , r.(i) 

/v£ — »v | T ^2 

z,(2) _ d 2 ) 4- t (2) 

h/r — /v | 1 n,2 

7 Numerical Examples 

Introduction 

The purpose of the numerical examples is: (i) to determine the influence of errors of 
alignment on the transmission errors and the shift of the bearing contact, and (ii) to prove 
that the predesigned parabolic function is able to absorb the transmission errors that are 
caused by the errors of alignment. We emphasize that the determination of transmission 
errors caused by misalignment is based on the following approach: 

(1) We consider an imaginary process for generation when an ideal transmission func- 
tion is provided. Then, using TCA, we simulate errors of alignment and determine the 
transmission errors that are caused by the respective error of alignment. 
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(2) On the second stage of investigation, we consider again the imaginary method of 
generation that has been provided in the previous section. We remind that this method 
of generation provides for each cycle of meshing a transmission function and a predesigned 
parabolic function as the sum of the ideal linear function and a predesigned parabolic func- 
tion. Using the TCA, we consider the meshing and contact of the gear misaligned gear drive 
that allows to determine: (i) the resulting function of transmission errors as the sum of the 
predesigned parabolic function and the function of errors that is caused by misalignment. 
We are also able to determine the shift of the bearing contact caused by misalignment by 
applying the TCA method. 

The simulation of meshing and contact has been accomplished for both methods of gen- 
eration described above that provide the longitudinal bearing contact, and the across the 
surface bearing contact. The results of computation confirmed that the bearing contact is 
stable, and the predesigned parabolic function is able indeed to absorb the almost linear 
functions of transmission errors caused by the respective errors of alignment. 

Input Data: 

The input data is represented in Tables 1-5. 

Output Data: 

The results of computation are represented for two cases of generation: (i) by application 
of a pinion head-cutter with straight blades (figs. 19-42); (ii) by application of a pinion 
head-cutter with circular arc blades (figs. 43-66). It is assumed that in both cases the gear 
is generated bv a head-cutter with straight blades. In each case, four sets of figures represent 
the respective influence of H , Q, V, S ' , which indicate the axial displacement of the pinion, 
the gear, the offset, and the change of the shaft angle, respectively. Alignment errors are 
given in millimeters, S' is given in arc minutes. 
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Table 1: Blank Data 



Pinion 

Gear 

Number of teeth 

11 

41 

Shaft angle 



Mean spiral angle 

35° 

35° 

Hand of spiral 

LF 

RH 

Outer cone distance (mm) 


90.07 

Face width (mm) 



Whole depth (mm) 

10.0 

10.0 

Pitch angle 


KeHaal 

Root angle 


74*59' 

Face angle 

15°T 

74°59' 


In each set of figures, we represent: 

(i) The initial influence of misalignment on the transmission errors, when the predesigned 
parabolic function has not been applied. The transmission function caused by misalignment 
is almost a linear function (see. for instance, figs. 19 and 20). 

(ii) The interaction of the linear function of transmission errors with the predesigned 
parabolic function. The results of TCA show that the obtained resulting function is indeed 
a parabolic function (see, for instance, figs. 21 and 23). 

(iii) The location of bearing contact for a misaligned gear drive (see, for instance, figs. 
23 and 24). 

Similar sets of figures are represented for other kinds of misalignment. 
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Table 2: Gear Cutter Specification 


Blade angle 

20° 

Cutter diameter (mm) 

152.4 

Point width (mm) 

2.79 


Table 3: Gear Machine Tool Settings 


Radial setting (mm) 

70.53744 

Cradle angle 

— 62‘ > 14 / 

Machine center to back (mm) 

0 

Sliding base (mm) 

0 

Blank offset (mm) 

0 

Machine root angle 

74°59' 


Table 4: Pinion Machine Tool Settings for Generation by a Cone 



Convex 

Concave 

Cutter blade angle 


EEHNil 

Cutter point radius (mm) 

71.7222 

80.4876 

Radial setting (mm) 

68.04991 

73.31925 

Cradle angle 

— 57 o 50' 

-66° 12' 

Machine center to back (mm) 

0 


Sliding base (mm) 

0 


Blank offset (mm) 

0 


Machine root angle 

15°!' 

i5°r 
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Table 5: Pinion Machine Tool Settings for Generation by Head-Cutter with circular Arc 
Blades 



Convex 

Concave 

Cutter point radius (mm) 

71.7222 

80.4876 

Radial setting (mm) 

68.04991 

73.31925 

Cradle angle 

—57*50' 

-66*12' 

Machine center to back (mm) 

0 

0 

Sliding base (mm) 

0 

0 

Blank offset (mm) 

0 

0 

Machine root angle 

15*1' 

15*1' 


8 Conclusion 

Extension of application of a CNC machine for generation of spiral bevel gears with the fol- 
lowing features has been discussed: (i) The gears are face-milled, the tooth depth is uniform, 
(ii) Two types of bearing contact are provided directed (a) in the longitudinal direction, and 
(b) in the direction across the surface, (iii) A predesigned parabolic function of transmission 
errors is provided for the absorption of transmission errors caused by misalignment, (iv) 
Equations of generated pinion-gear tooth surfaces have been derived, (v) TCA computer 
program has been developed and the influence of misalignment on the transmission errors 
and the shift of the bearing contact has been investigated. The computations that were 
performed confirmed the stability of the bearing contact, the low level of transmission errors 
and the favorable shape of the function of transmission errors, of a parabolic type, for a 
misaligned gear drive. 
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9 Directions for TCA Program Use 

There are two TCA programs for two different cases. One program is BEVEL.FOR, for the 
case that both the pinion and the gear generating surfaces are cone surfaces. The other 
program is RBEVEL.FOR, for the case that the pinion generating surface is a surface of 
revolution, the gear generating surface is a cone surface. For both programs the input data 
files and the output data files are almost the same, except that there is an additional arc 
radius, RHO, in the input data files for program RBEVEL.FOR. 

Input data 

1. Control codes 

(a) For right hand gear JCH=1, for left hand gear JCH=2 

(b) TL1 and TL2 are numbers of extra point on the contact path which should not be 
larger that 2 

(c) MM is the number of contact points 

2. Blank data 

TNI — Pinion number of teeth 

TN2 — Gear number of teeth 

C — Shaft offset (zero for spiral bevel gear) (mm) 

TW — Face width of gear (mm) 

GAMMA — Shaft angle (degree) 
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MCD — Mean cone distance (mm) 

RGMA1 — Pinion root cone angle (degree) 

B1 — Pinion spiral angle (degree) 

B2 — Gear spiral angle (degree) 

RGMA2 — Gear root cone angle (degree) 
FGMA2 — Gear face cone angle (degree) 
PGMA2 — Gear pitch cone angle (degree) 
ADD2 — Gear mean addendum (mm) 

DED2 — Gear mean dedendum (mm) 

WD — Whole depth (mm) 

CC — Clearance (mm) 

DEL — Elastic approach (mm) 

3. Gear cutter specification 

RU2 — Gear nominal cutter radius (mm) 
PW2 — Point width of gear cutter (mm) 
ALP2 — Blade angle of gear cutter (degree) 

4. Gear machine-tool settings 

XG2 — Machine center to back (mm) 
GAMA2 — Gear machine root angle (degree) 
XB2 — Sliding base (mm) 
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EM2 — Blank offset (mm) 


4. Pinion machine-tool settings 

RCF — Point radius (mm) 

XG1 — Machine center to back (mm) 

XB1 — Sliding base (mm) 

EMI — Blank offset (mm) 

GAMA1 — Pinion machine root angle (degree) 

A1P1 — Blade angle of pinion cutter (degree) 

6. Misalignments 

A — Constant coefficient of the predesigned parabolic function 
H — Misalignment along the pinion axis (mm) 

Q — Misalignment along the gear axis (mm) 

V — Misalignment of axis offset (mm) 

S ' — Misalignment of shaft angle (arc min.) 

Input data files 

Files 70 and 60 are for program RBEVEL.FOR, file 70 for the convex side, file 60 for 
the concave side. These two files must be read together. Files 90 and 80 are for program 
BEVEL. FOR, file 90 for the convex side, file 80 for the concave side. These two files must 
be read together. 

Output data files 
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File 9 is an overall output data file. All the input and output information is stored in this 
file. File 91 stores the information of transmission errors for the convex side; file 93 stores 
the information of transmission errors for the concave side. File 92 stores the information of 
contact path and contact ellipse for the convex side; file 94 stores the information of contact 
path and contact ellipse for the concave side. 
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Figure 3: Generating surfaces E c and E ( that provide contact path along the surface 
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Figure 4: Generating surfaces £ c and £ t that provide a contact path across the surface 


41 



Figure 5: Generating blades and cones 
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Figure 8: Machine-tool settings for gear generation 
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Figure 9: Gear generating surface and its straight line blade profiles 
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Figure 12: For derivation of relation between imaginary and real pinion machine-tool settings 
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Convex 


Figure 13: Profiles of pinion and gear head-cutters represented in plane Z m - 0 
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( a) For concave side 
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(b) For convex side 


Figure 14: Profiles of straight-line pinion blades 
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Figure 15: For derivation of pinion generating cone surface 
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<t > i (deg.) 


Figure 20: Transmission errors: H=0.01mm (Concave side) 
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Figure 21: Resulting transmission errors: H=0.01mm (Convex side) 



Figure 22: Resulting transmission errors: H=0.01mm (Concave side) 
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Tooth length (mm) 


Figure 23: Shift of bearing contact: H=0.01mm (Convex side) 



Tooth length (mm) 


Figure 24: Shift of bearing contact: H=0.01mm (Concave side) 
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Figure 27: Fesulting transmission errors: Q=0.01mm (Convex side) 



Figure 28: Resulting transmission errors: Q=0.01mm (Concave side) 
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Tooth height (mm) 



Tooth length (mm) 


Figure 29: Shift of bearing contact: Q=0.01mm (Convex side) 



Tooth length (mm) 


Figure 30: Shift of bearing contact: Q=0.01mm (Concave side) 
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Figure 31: Transmission errors: V=0.01mm (Convex side) 



Figure 32: Transmission errors: V=0.01mm (Concave side) 
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4 > i (deg.) 


Figure 33: Resulting transmission errors: Y=0.01mm (Convex side) 



4 > i (deg.) 
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Figure 34: Resulting transmission errors: V— 0.01mm (Concave side) 



Figure 35: Shift of bearing contact: V=0.01nun (Convex side) 



Figure 36: Shift of bearing contact: V=0.01mm (Concave side) 
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Figure 38: Transmission errors: <5'=1 arc min. (Concave side) 
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Figure 41: Shift of bearing contact: £'=1 arc min. (Convex side) 



Tooth length (mm) 


Figure 42: Shift of bearing contact: <5'=1 arc min. (Concave side) 
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Fignre 43: Transmission errors: H=0.01mm (Convex side) 
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Figure 44: Transmission errors: H=0.01mm (Concave side) 
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4>i (deg.) 


Figure 45: R°sulting transmission errors: H=0.01mm (Convex side) 




Figure 47: Shift of bearing contact: H=0.01mm (Convex side) 



Tooth length (mm) 


Figure 48: Shift of bearing contact: H=0.01mm (Concave side) 


70 



Figure 49: Transmission errors: Q=0.01mm (Convex side) 



Figure 50: Transmission errors: Q=0.01mm (Concave side) 
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Figure 51: Resulting transmission errors: Q=0.01mm (Convex side) 



Figure 52: Resulting transmission errors: Q=0.01mm (Concave side) 
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Figure 53: Shift of bearing contact: Q=0.01mm (Convex side) 
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Figure 54: Shift of bearing contact: Q=0.01mm (Concave side) 
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Figure 57: Resulting transmission errors: V=0.01mm (Convex side) 



Figure 58: Resulting transmission errors: V=0.01mm (Concave side) 
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Figure 59: Shift of bearing contact: V=0.01mm (Convex side) 



Tooth length (mm) 


Figure 60: Shift of bearing contact: V=0.01mm (Concave side) 
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Figure 64: Resulti 



Figure 65: Shift of bearing contact: 6 '= 1 arc min. (Convex side) 



Figure 66: Shift of bearing contact: £'=1 arc min. (Concave side) 
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